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The ground-state properties of the S — 1/2 frustrated Heisenberg spin chain with interactions 
up to fourth nearest neighbors are investigated by the exact-diagonahzation method and density 
matrix renormahzation group method. Our numerical calculations clarify that the ferrimagnetic 
state is realized in the ground state in spite of the fact that a multi-sublattice structure in the 
shape of the system is absent. We find that there are two types of ferrimagnetic phases: one 
is the well-known ferrimagnetic phase of the Lieb-Mattis type and the other is the nontrivial 
ferrimagnetic phase that is different from that of the Lieb-Mattis type. Our results suggest 
that a multi-sublattice structure of the shape is not necessarily required for the occurrence of 
ferrimagnetism. 
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Ferrimagnetism is one of fundamental phenomena in 
the field of magnetism. A typical case showing ferrimag- 
netism is that when a system includes spins of two types 
that antiferromagnetically interact between two spins 
of different types in each neighboring pair. The sim- 
plest example is an (5, s)=(l, 1/2) antiferromagnetic 
mixed spin chain, in which two different spins are ar- 
ranged alternately in a line and coupled by the nearest- 
neighbor antiferromagnetic interaction.-'^' The occurrence 
of ferrimagnetism in this case is understood within the 
Marshall-Lieb-Mattis theorem concerning quantum spin 
systems.^' '^^ Even though a system includes spins of one 
type, this theorem also derives the presence of ferrimag- 
netism when the system includes more than one sub- 
lattice of spin sites, for example, the spin system in a 
diamond chain. From these two mechanisms, the ex- 
istence of a multi-sublattice structure is very important 
for the occurrence of ferrimagnetism. 

At this stage, one asks a fundamental question: Is a 
multi-sublattice structure in the shape of a Hamiltonian 
essential and necessary for the occurrence of ferrimag- 
netism? The purpose of the present study is to answer 
this question. Our following demonstration will clarify 
that the answer is no. In this study, wc find that ferri- 
magnetism can appear due to the effect of magnetic frus- 
tration even in the absence of a multi-sublattice structure 
in the shape of a system. 

In this study, we examine the model whose Hamilto- 
nian is given by 



n 



1 



(1) 



Si-i-3 + 2^^* 



Si+2 + Si • Si+4)], 



where Si is the 5* = 1/2 spin operator at the site i. The 
system size is denoted by N. We emphasize here that this 
model has only one spin in a unit cell, namely, it has no 
sublattice structure. Energies are measured in units of 
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J; therefore, we set J = 1 hereafter. We have a control- 
lable parameter, J', in the Hamiltonian (1). This model 
was originally introduced in ref. 11 detailing the study 
of constructing a model Hamiltonian as a generalization 
from the Majumdar- Ghosh model. The Hamiltonian 
(1) includes two cases in which the ground state of the 
system is exactly obtained. For J' — 0, the system is re- 
duced to the Majumdar- Ghosh model, whose ground 
state is described by direct products of spin-singlet states 
in nearest-neighbor pairs of S" = 1/2 spins. The ground 
state is called the dimer (DM) state. Note that even if 
J' takes a nonzero value, this DM state is still an eigen- 
state of the system. The DM state becomes an excited 
state when J' increases. In the limit of a large J', on 
the other hand, the ferromagnetic (FM) state becomes 
the ground state. Although the wavefunctions of these 
limits are well known, the ground state in the interme- 
diate region is not sufficiently understood. In ref. 11, it 
was reported that the spontaneous magnetization in the 
intermediate region appears and that the magnetization 
changes gradually. In the present study, we investigate 
the magnetic structure of the ground state in this inter- 
mediate region by some numerical calculations. We show 
that our results lead to the conclusion that the ferrimag- 
netic state can appear in the ground state, even of models 
consisting of only a spin in each unit cell. 

We employ two reliable numerical methods, the ex- 
act diagonalization (ED) method and density matrix 
renormahzation group (DMRG) method."'") ^he ED 
method can be used to obtain precise physical quanti- 
ties for finite-size clusters. This method does not suffer 
from the limitation of the shape of the clusters. It is 
applicable even to systems with frustration, in contrast 
to the quantum Monte Carlo (QMC) method coming 
across the so-called negative-sign problem for a system 
with frustration. The disadvantage of the ED method 
is the limitation that the available sizes are only small. 
Thus, we should pay careful attention to finite-size ef- 
fects in quantities obtained from this method. On the 
other hand, the DMRG method is very powerful when a 
system is one-dimensional under the open-boundary con- 
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Fig. 2. (Color) Size dependences of the phase boundaries. The 
results presented are those oi N = 12, 18, 24, and 30 from the ED 
calculations and those of N = 48, 72, and 96 from the DMRG 
calculations. Squares (circles) denote results in the cases under 
the periodic (open)-boundary condition. Dotted lines are drawn 
as guides for the eyes between the data from the ED and DMRG 
calculations. In the limit N — > oo, the phase boundary between 
the < M/Ms < 1/3 phase and the M/Ms = 1/3 phase seems 
to converge to approximately J' = 1.30. 




Fig. 1. (Color) (a) Lowest energy in each subspace divided by 
51°*. Results of the DMRG calculations are presented when the 
size system is Af = 96 for J' = 0.685,0.87, and 1.5. Arrows in- 
dicate the spontaneous magnetization M for a given J'; M is 
determined to be the highest 5*°* among the values taking the 
lowest common energy, (b) J' dependence of the normalized mag- 
netization M/Ms in the ground state. Red circles (black squares) 
denote the results obtained by ED calculations for a size system 
of Af = 24 under the open (periodic)-boundary condition. In the 
inset of (b), blue diamonds show the results obtained by DMRG 
calculations for a size system of A'^ = 96 under the open-boundary 
condition accompanied by red circles denoting the results ob- 
tained by ED calculations for a size system of TV = 24 under the 
open-boundary condition. 

dition. The method can treat much larger systems than 
the ED method and is apphcable even to a frustrated sys- 
tem. In the present research, we use the "finite-system" 
DMRG method. 

In the present study, two quantities are calculated by 
the two methods mentioned above. One is the lowest en- 
ergy in each subspace divided by 5*°' to determine the 
spontaneous magnetization M, where 5'*°' is the z com- 
ponent of the total spin. We can obtain the lowest energy 
E{N, Sl°^, J') for a system size N and a given J'. For ex- 
ample, the energies of each 5*°* in the three cases of J' 
are presented in Fig. 1(a). This figure is obtained by our 
DMRG calculations of the system of iV = 96 with the 
maximum number of retained states {MS) 1500, and a 
number of sweeps [SW] 50. The spontaneous magneti- 
zation M for a given J' is determined as the highest 5"*°' 
among those at the lowest common energy. (See arrows in 
Fig. 1(a).) The other quantity is the local magnetization 
in the ground state for investigating the spin structure 
of the state. The local magnetization is obtained by cal- 
culating {Sf), where {A) denotes the expectation value 
of the physical quantity A and S^ is the z-component of 
the spin at the site i. 

First, let us examine the J' dependence of M/M^ to 
confirm the existence of the intermediate phase between 
the FM phase and the nonmagnetic DM phase irrespec- 
tive of the boundary conditions, where Mg is the satu- 
ration value of the magnetization. Results are presented 



for iV = 24 from our ED calculations under the open 
and periodic boundary conditions in Fig. 1(b). We suc- 
cessfully observe the intermediate-magnetization phase 
irrespective of the boundary conditions. We also include 
in Fig. 1(b) some DMRG results of TV = 96, which sug- 
gests a weak size dependence of M/Mg as a function of 
J'. Careful observation of the region of < M/Ms < 1/3 
enables us to find that the intermediate-magnetization 
phase consists of two phases. One is the phase where 
M/Ms is fixed at 1/3; this feature is that of the ferri- 
magnetism of the so-called Lieb-Mattis (LM) type, in 
which the spontaneous magnetization is fixed to be a 
simple fraction of the saturated magnetization.^''^) The 
other is the phase where M/Ms changes continuously 
with respect to the strength of J'. This feature is cer- 
tainly different from that of the LM ferrimagnetism; the 
continuous change in M/Ms is observed as the ferrimag- 
netism of the non-Lieb-Mattis (NLM) type in several 
models. ^^"^^) We will determine later whether or not the 
phase of < M/Ms < 1/3 in the present model is of the 
NLM type. Note here that these two phases are observed 
under both boundary conditions. On the other hand, the 
region of 1/3 < M/Ms < 1 is observed near M/Ms = 1 
only under the open-boundary condition. At present, it 
is unclear whether or not this phase survives in the limit 
oo. 

Next, we study the size dependences of the bound- 
aries between the phases observed above. We investigate 
five boundaries: J' = J( between the DM phase and the 
phase of < M/Ms < 1/3, J' = J2 between the phase 
of < M/Ms < 1/3 and the phase of M/Ms = 1/3, 
J' = J3 between the phase of M/Ms = 1/3 and the 
phase of 1/3 < M/Ms < 1, J' — J4 between the phase 
of 1/3 < M/Ms < 1 and the FM phase, and J' = J'^ 
between the phase of M/Ms = 1/3 and the FM phase 
without the phase of 1/3 < M/Ms < 1- Note that J3 and 
J4 appear under the open-boundary condition, whereas 
Jg appears under the periodic-boundary condition. Fig- 
ure 2 shows the results of = 12, 18, 24, and 30 from 
the ED calculations and those of A^ = 48, 72, and 96 
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Fig. 3. (Color) Local magnetization (5"^^) under the open- 
boundary condition: (a) for J' = 1.5, (b) for J' = 0.87, and 
(c) for J' = 0.685 from the DMRG calculation for N = 96. The 
site number is denoted by i, which is classified into i = 3ra — 2, 
3ti — 1, and 3n, where n is an integer. Squares, circles, and tri- 
angles mean i = 3n — 2, 3n — 1, and 3n, respectively. 



from the DMRG calculations. One finds that J[ from the 
ED calculations under the periodic-boundary condition 
and that from the DMRG calculations under the open- 
boundary condition are consistent with each other; we 
have J[ ^ 0.59 as an extrapolated value. Concerning the 
boundary Jj, there exists a not so small difference be- 
tween the result under the open-boundary condition and 
that under the periodic-boundary condition for a given 
N; however, seems to converge to 1.30 irrespective of 
the boundary condition. On the other hand, the situa- 
tions of the boundaries of the phase of M /Ms = 1/3 and 
the FM phase are slightly complicated in our results. It 
seems that Jg and J4 become farther away from each 
other with increasing N and that J3 and J5 converge to 
the same value of 1.77. We also have J4 converging to 
2.06. From these results of the extrapolation, it is evi- 
dent that the phase of M/Mg = 1/3 and the phase of 
< M/Ms < 1/3 exist in the thermodynamic limit. On 
the other hand, it is difficult to determine whether or 
not the phase of 1/3 < M/Mg < 1 is present. There is a 
possibility that this phase merges with the FM phase in 
the thermodynamic limit for two reasons: one is that this 
phase appears only near M/M^ — 1 and the other is that 
it is observed only under the open-boundary condition. 
The issue of whether or not this phase survives should 
be clarified in future studies; hereafter, we do not pay 
further attention to this phase. 

Next, we examine the local magnetization (5f ) in the 
two phases of < M/M^ < 1/3 and M/M^ = 1/3 to 
determine the magnetic properties in each phase. We 
present our DMRG results of (S'f) of the system of 




Fig. 4. (a) Spin configuration from the point of view of classical 
vectors. The site number i in the Hamiltonian (1) is classified 
into 3n, 3n — 1, and 3n — 2, where n is a positive integer. The 
angle 6 for J' is determined by minimizing the classical energy, 
(b) J' dependences of classical energies of eq. (2) for 8 = (LM, 
dotted line), it/3 (NM, dotted chain line), and FM (solid Une) 
energy of eq. (3). 



= 96. Note here that we calculate (S^) within the 
subspace of the highest 5*°* corresponding to the spon- 
taneous magnetization M obtained for a given J'. The 
results of (S'f ) are shown in Figs. 3(a)-3(c) for J' = 0.685, 
0.87, and 1.5, respectively. In each case, one can observe 
a three-sublattice structure of the spin state clearly. In 
Fig. 3(a), the i dependence of (S'f ) in each of the sub- 
lattices of the spin structure is weak around the center 
of the system, although the edge effect spreads into a 
wide range from the edges. This behavior suggests that 
the spin state forms the LM ferrimagnetic state of up- 
up-down, which is consistent with AI/Mg—l/S in the pa- 
rameter region near approximately J' — 1.5. In Figs. 3(b) 
and 3(c), on the other hand, we find that the local mag- 
netization shows a longer-distance periodicity in addition 
to the three-sublattice structure. The longer-distance pe- 
riodicity changes when J' is changed within the phase of 
< M/AIs < 1/3, the periodicity suggests an incom- 
mensurate modulation. A similar feature of this local 
structure was reported in some one-dimensional quan- 
tum frustrated spin systems. -'^^^ Therefore, the phase 
of < M/Ms < 1/3 is considered as the NLM-type ferri- 
magnetic phase. This incommensurate feature originates 
from the effects of quantum fluctuation and frustration. 
We also calculate {S^) for different system sizes, N = 48 
and 72. At least from these data (not shown in this pa- 
par), the periodicity and amplitude of the modulation 
seem to show only weak dependences on the system size. 
Note that the behavior of long-distance periodicity ac- 
companied by the three-sublattice structure at the same 
time is different from the wave functions with a long pe- 
riodicity reported in ref. 23. 

Here, let us discuss the behavior of the intermediate 
phase between the FM phase and the nonmagnetic phase 
from the viewpoint that spins in the Hamiltonian (1) are 
assumed to be classical vectors. We consider the spin con- 
figuration of the classical vectors depicted in Fig. 4(a), 
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where the characteristic angle is defined. This classical 
spin arrangement has been determined from our obser- 
vation in Fig. 3 that the three-sublattice spin structure 
is realized in the intermediate region. The case oi 9 = 
means that this classical state is the LM-type ferrimag- 
netic state with the ratio of the spontaneous magnetiza- 
tion to the saturated magnetization to be 1/3. On the 
other hand, 9 = ir/'i means that the state is in a non- 
magnetic (NM) state. The classical energy per spin site 
under the periodic-boundary condition is given by 

E{J',e) = ^[{&-AJ')cos^e-{&-AJ')cos9+{-2,-AJ% 

(2) 

and the energy of the ferromagnetic state is given by 



FM 



(3-4J')/8. 



(3) 



The dependences of the energies shown in eqs. (2) and 
(3) are shown in Fig. 4(b). The FM (NM) phase appears 
at J' > 1.5 (J' < 1.5). One finds that J'=1.5 is the 
boundary of the FM and NM phases. At exactly J' =1.5, 
many states degenerate, including not only the FM and 
NM states but also the fcrrimagnetic state with an ar- 
bitrary angle 9. There is no intermediate phase between 
the two phases. It is worth emphasizing here that even 
the LM fcrrimagnetic phase does not appear. This argue- 
ment suggests that the occurrence of the intermediate- 
magnetization state observed in the Hamiltonian (1) of 
the quantum system is a consequence of the quantum 
effect induced by frustration. 

Finally, we mention another case when the 
intermediate-magnetization phase appears in the 
frustrated spin system in one dimension with anisotropic 
interactions. Note here that this phase disappears 
in the isotropic case of interactions, which suggests that 
the origin of this phase is the anisotropy. However, it 
has not been examined yet whether or not this model 
shows a similar incommensurate modulation. Such 
examination would clarify the relationship between the 
intermediate magnetization of this model and the NLM 
ferrimagnetism studied in the present case. 

In summary, we study the ground-state properties 
of an S = 1/2 frustrated Heisenberg spin chain with 
isotropic interactions up to the fourth nearest neighbor 
by the ED and DMRG methods. In spite of the fact that 
this system consists of only a single spin site in each unit 
cell determined from the shape of the Hamiltonian, the 
ferrimagnetic ground state is surprisingly realized in a fi- 
nite region between the ferromagnetic and nonmagnetic 
states. This result is in contrast to that of other systems 
of translationally invariant chains. We find that the 
intermediate region consists of phases of two ferrimag- 
netic types, the Lieb-Mattis type and non-Licb-Mattis 
type. In the latter phase, we confirm that the local mag- 
netization shows characteristic incommensurate modula- 
tion. The presence of the fcrrimagnetic state without a 
sublattice structure of the shape of the system is a con- 
sequence of the strong quantum effect induced by frus- 
tration. Our findings shed light on a new aspect of the 
effect of frustration in quantum systems. 
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